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Abstract 

We study classical energy-momentum (E-m) diagrams for rotational motion of dipolar asymmet- 
ric top molecules in strong external fields. Static electric fields, nonresonant linearly polarized laser 
fields, and collinear combinations of the two are investigated. We treat specifically the molecules 
iodobenzene (a nearly prolate asymmetric top), pyridazine (nearly oblate asymmetric top), and 
iodopentafluorobenzene (intermediate case). The location of relative equilibria in the E-m plane 
and associated bifurcations are determined by straightforward calculation, with analytical results 
given where possible. In cases where analytical solutions cannot be obtained, we resort to numer- 
ical solutions, while keeping a geometrical picture of the nature of the solutions to the fore. The 
classification we obtain of the topology of classically allowed rotor configuration space regions in 
the E-m diagram is of potential use in characterization of energy eigenstates of the corresponding 
quantum mechanical problem. 
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1. INTRODUCTION 



The study of rigid body motion is one of the most important topics in classical and 
quantum mechanics l|, 1^, 4, SjIgI Q, Q]. Of particular importance are the integrable 



which include the free asymmetric top (Euler top ) 
4, 6], the symmetric top in a uniform external gravitational field (Lagrange top) [g], 



cases of the rigid body problem [9 
1 

and the Kovalevskaya top [6 



12l | . The theory of rigid body motion also provides the basis for 



analysis and interpretation of the rotational dynamics and spectra of semi-rigid molecules 



0,y 



13 



molecule 



14 



14 



161 ] ■ In molecular terms the Euler top is simply a free asymmetric top 



151 ]. the Lag r ang e top models a symmetric top molecule with a dipole 



moment in an electric field 



15 



17| , while there does not appear to be an obvious molecular 



analogue for the Kovalevskaya top. Breaking the symmetry of the moment of inertia tensor 
in the Lagrange top results in the non-integrable problem of an asymmetric to p in a static 
field. Molecular examples of these two cases are iodobenzene (near prolate) 18|, ll9| and 



pyridazine (near oblate) j20|] in a static electric field. The rotational constants of the molecule 



iodopentafiuorobenzene make it a more generic example of an 



asymmetric rotor 



2l|. 



The general problem of classical-quantum correspondence |22l . 1231 ] is of great interest for 
both integrable and nonintegrable rotor systems. Various aspects of the classical-quantum 
correspondence have been studied for diatomic molecules in tilted fields, i.e., noncoUinear 
static electric and nonresonant linearly polarized laser fields 2J]. The integrable coUinear 
case exhibits the phenomenon of monodromy 11[ both classically and quantum mechanically 



251 ]. For the nonintegrable case of tilted fields the rotor motion tends to be integrable in 



both low-energy (pendular) and high-energy (free-rotor) limits, and chaotic at intermediate 
energies, with the degree of chaos controllable by variation of the angle between the fields 



241 ]. For collinear fields the system is integrable, with both the energy E and the projection 



of the angular momentum into the space fixed z-axis, m, as constants of motion. The effective 



potential Vcs{0] rn) for a given value o: 
define the relative equilibria 



' m exhibits extrema in the 6 (polar) coordinate, which 
261 1. Plotting the location of these extrema in the E-m 

n 

plane gives the energy- momentum diagram for the system [27[ . The E-m diagram provides 
a useful global classification of the rotor dynamics, as distinct regions of the E-m plane are 
associated with different allowed types of motion of the diatomic. For symmetric tops in 
electric fields similar diagrams can be constructed [l3], and analysis of classical symmetric 



2 



;op E-m and E-k diagrams helps understand the organization of the quantum level spectrum 



The problem of a dipolar asymmetric top in a static external field is nonintegrabl e |27 .l28| 



as is the problem of a polarizable asymmetric top in a nonresonant laser field [27|. For 
asymmetric tops in either static or laser fields or collinear superpositions of the two, the 
angular momentum projection m is a constant of motion. Although the complicated form of 
the kinetic energy does not allow us to separate an effective potential as straightforwardly as 
in the diatomic or the symmetric top case, it is still possible to define an effective or amended 



potential [6|, |10| 



2^ for the class of motions in which the asymmetric top molecule rotates 
with constant Euler angles 6 and ip (the third Euler angle is an ignorable coordinate). 
The energies of the extrema of this potential for given values of m again define an energy- 
momentum diagram ^, 26, 29, 30, 31; 32|, which can be used to classify the motions of the 
asymmetric top. 

The asymmetric top in an external static electric field is an example of a dynamical 



system with symmetry 



33, 



3J] . The potential energy and Hamiltonian in this case are 



invariant with respect to rotations about the space-fixed field direction. There is therefore an 
associated constant of the motion m, the projection of the angular momentum vector j onto 
the external field axis, in addition to the energy E. The mapping of the system phase space 



onto the E-m plane is called the energy-momentum map [61, 1261 



33 



341 ]. Critical points of this 



mapping define the bifurcations sets in the E-m plane, which form the boundaries of regions 
of qualitatively different types of classical motion l3, [26]. Following the fundamental work 
of Smale fl. appl.cations of these coneepts were made to integxable rotor problems fl, 
asymmetric rotors in an external gravitational field [30| , and symmetric [35| and asymmetric 



3l| rotors in more complicated potentials. Relative equilibria have a 



rotating semi-rigid molecules in the absence of external fields 



36 



so been studied for 



states in rotationally inelastic collisions 38 



37( 1 and for transition 



In the present paper we study E-m diagrams of relative equilibria for molecular asymmet- 
ric tops for several molecule- field configurations of physical interest 39|, |40f|: a static electric 
field; a nonresonant linearly polarized laser field; both fields in a collinear combination. For 
the field strengths considered, the free asymmetric top motion is strongly perturbed. The 
associated E-m diagrams are obtained for the most part analytically. Our aim here is to 
understand the nature of the rotor motions associated with different regions and curves in 
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the E-m diagrams for physically relevant values of the field parameters. We also study the 
topological classification of the allowed O-ip configuration space of the system in terms of 
their (multivalued) genus [6|. An extension of these results, to be discussed in a future paper, 
involves comparison of computed quantum mechanical eigenstate probability densities with 
the boundaries of classically allowed regions in O-ip configuration space. The comparisons 
indicate that the classical mechanical methods developed here provide a promising founda- 
tion for the difficult task of classifying the quantum levels of the complex system consisting 



of an asymmetric rotor in external fields 



41 



42, 



43|. 



We mention that the approach adopted here may be thought of as the analogue for 



perturbed rotor systems of the methods applied by Kellman and coworkers 4^, |45| to vibra- 



tional problems exhibiting a single conserved vibrational (superpolyad) quantum number 



fsee also Ref. 



46|). 



There have been many studies of the classical dynamics of a rigid asymmetric top rotating 



about a fixed point in a gravitational 



example 



m 



30 



32 



471, 



48 



49 



fie 



50 



d (this is the classic heavy top problem); see, for 



5l| . Katok [30|] and Gashenenko and Richter 51 1 



obtained E-m diagrams and analyzed bifurcations of relative equilibria in the E-m plane. 
These authors also classified the topology of accessible system configuration space. We 
present here a similar analysis for several molecular examples of an asymmetric top molecule 
possessing a dipole moment in a static electric field. In the spirit of Katok's analysis jsol , 
we use straightforward analytical and geometric methods to build the E-m diagram for the 
molecules of interest. 

Katok's treatment was generalized by Tatarinov 311] (see also 6|]) to include more compli- 
cated gravitational perturbations of the rotational dynamics of an asymmetric top. Although 
given for a specific potential, Tatarinov's analysis can be mapped directly onto the molecu- 
lar case of an asymmetric top in collinear fields 27|]. Following Tatarinov, we obtain E-m 
diagrams for asymmetric molecules of physical interest, and study bifurcations and other 
characteristic features of the problem, including the topology of the 6-ip configuration space 
and its classification according to the genera of the (connected) allowed regions Q]. 

This paper is organized as follows. First in Section [2] we derive the Hamiltonian for 
asymmetric tops in the general tilted fields case, with the aim of introducing our notation 
and conventions. In Section [3] we treat the molecules in a static electric field, while in Section 
m we study the same molecules in a nonresonant linearly polarized laser field. In Section [5] 



4 



we analyse these molecules in coUinear fields. Section [6] concludes. 

Finally we mention that, the advantages of alternative approaches notwithstanding 52| . 
all calculations reported here have been carried out using polar coordinates. 
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2. HAMILTONIAN FOR THE ASYMMETRIC TOP IN COMBINED FIELDS 



2.1. General case: tilted fields 

We use the ?/-convention to define the three Euler angles {6, (j), ifj) describing the 
orientation of the body-fixed frame with respect to lab-fixed frame. In the body-fixed frame 
the kinetic energy of the free asymmetric top can be written in terms of the components of 
the angular momentum j = (ji, j2, Ja) and the three components of the (diagonal) moment 
of inertia tensor I = diag(/i, I2, 13) 



T 



2/1 



(2.1) 



2/2 2J3 

For an asymmetric top, Ii ^ I2 ^ Is- In terms of the Euler angles and their conjugate 
momenta {pe,P(t>,Pij) the body-fixed components of j are 

cosip 



Ji 

32 



sin — J 

pe cos ip + 
P^ 



' sin 9 

sm.1l) 

^ sin 6 



+ p^ cos ip cot 6*, 
— p^ sin ip cot 6, 



(2.2a) 
(2.2b) 
(2.2c) 



so that 



■2 -2 , -2 , -2 

J = Jl + J2 + Js 



sin'^ 



p^ cos9f +pI. 



From (12. ip . the kinetic energy T = T{6, </>, ip,pe,P(i>,P^) is then 

2 



T 



2/1 



COS Ip 

Pe sin if) + {p^ cos 6 - p^) 



1 

2^2 



Pe cos i) + ^— 
sm 6 



p^ cos I 



n 2 



(2.3) 



j4 
2/3' 

(2.4) 



Immediately we see that is an ignorable coordinate and p^ = m, the projection of j into 
the space fixed z-axis, is a constant of the motion for the free top. 

If the polarizability in the molecule-fixed frame is given by the diagonal tensor a = 



diag(ai, 02,03! 
fixed 2;-axis is 



jthe interaction with a nonresonant laser field polarized along the space- 



39| 



Vl = — ^ [ai + (0^2 — «i) sin^ 6 sin^ tp + {a^ — ai) cos^ 6*] 



(2.5) 
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with ef^ proportional to the intensity of the laser field. Omitting the angle-independent term 
(which produces a constant shift in energy) we obtain 

Vl = —Auj2 sin 6 sin V; - Acu3 cos^ ^, (2.6) 

with Au!2 = («2 — tti) £^i/4 and Au^ = (as — ai) Assuming the dipole moment to lie 

along the molecule-fixed z-axis, the interaction with a static electric field tilted through an 
angle f3 with respect to the space-fixed 2;-axis and lying in the space-fixed xy-plane is 

Vs = —doEsicos P cos 9 + sinP cos sin 9), (2.7) 

where do is the magnitude of the electric dipole moment and £5 is the strength of the static 
field. 

The Hamiltonian of the asymmetric top in tilted fields can be written from equations 
(El, (EZD and (El]) as 

H = T + Vs + Vl. (2.8) 

Note that the potential V5 -|- is a function of the angle 0, so that m is not conserved. The 
asymmetric top in tilted fields is therefore a physically significant rotor problem with three 
degrees of freedom. For collinear fields, /9 = 0, the angle is not present in H and Prp = m 
is then a constant of the motion. In the remainder of this paper we consider the collinear 
case only. 



2.2. Collinear fields 

The Hamiltonian for an asymmetric top molecule in collinear fields can be rewritten as 

H = li\-'j)-j + Vi9,^), (2.9) 

with 

V{9, ip) = —uj cos 9 — Auj2 sin 9 sin - Auo^cos^ 9. (2.10) 

In terms of the Euler angles and their time derivatives the body-fixed components of the 
angular momentum are 

ji = Ii (j)smip — (f> sin 9 cos , (2.11a) 
j2 = h (9cosip + 4>sm9smip^ , (2.11b) 
h = h(<Pcos9 + ij) . (2.11c) 
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The projection of j onto the space-fixed unit vector = (0,0, 1) is given in terms of the 
direction cosine matrix C as 

m = C3-j, (2.12) 

where 

C3 = (C31, C32, C33) 

(2.13) 

= (— sin ^ cos?/;, sin 6' sin?/;, cos 6'), 
is the vector of body-fixed components of e^. 



2.3. Determination of relative equilibria 



As the angle is ignorable, we consider relative equilibria defined by the conditions 6 = 0, 



and = 



30|]. These relative equilibria in general define periodic orbits in the full rotor 



phase space. Using these conditions in equations (12. lip , and rewriting in terms of G3 we 
obtain 

j = 0(IC3). (2.14) 



The equation for m, (12.121) . can be rewritten as 

m = 0(103) ■ Cs 



(2.15) 



which can be used to express in terms of m and C3. Substituting the resulting expression 
for j, 



■ _ m{\Cs) 

^-(1C3)-C3' 

into the Hamiltonian (12.91) gives the effective or amended potential [6|, llOl . 
asymmetric top molecule in coUinear fields 

,2 



(2.16) 



30 



3l[ | for an 



m 



(2.17) 



2(IC3) ■ C3 

Relative equilibria with 6 = 0, ip = are found from the effective potential (I2.17P solving 
the equations 



de 

dVm 

dijj 



[(IC3) ■ C;]^ de de 



0, 



[(IC3) • C3]2 di; 



^ = 0, 



(2.18a) 
(2.18b) 
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with 



dC 

-—— = (— cos ip cos 6', sin tp cos 6, — sin 9), (2.19a) 



dip 



(sin?/'sin6', cos?/'sin6', 0) (2.19b) 



and 



dV 

— — = cij sin 6' — 2 Au;2 sin 6 cos 6 sin^ ih + 2Au;3 sin 6 cos 6', (2.20a) 
dV 

— = -2^uj2siY).%ljcos%IJsm^ e. (2.20b) 

In general the solutions of equation fl2.18p depend on m. For a given m the zero set of 
(respectively, ^^) in the O-ip space gives the solution set for equation fl2.18al) (respectively, 
(]2.18b[) ). The solution set typically consists of (possibly disjoint) curves in the 6-ip plane. 
In practice (see below), we are able to give a natural parametrization of each of these 
solution curves, so that the general solution of equations (12.181) is obtained by evaluating 
(respectively ^^) along each of the particular solution curves. Zeroes of the relevant 
functions are found by interpolation along the curve, yielding values of 9 and ip that solve 
()2.18p for given m. Finally, the effective potential (12.171) is evaluated at each of these solution 
points, which in general have different energies for a given m. The calculation is repeated 
for different values of m to obtain the complete E-m diagram. 

To clarify the general procedure just outlined, consider the Euler problem (free asym- 
metric top), with u) = Auj2 = Auj^ = 0. Equations (I2.18P in this case are 



dVm sin 9 cos 9 (Ji cos^ i' + h sin^ ip — I3) 



m 



2 



89 [(IC3) ■ C3]2 

dVjn simp cosip sin^ 9 {I2 — Ii)m'^ 



0, (2.21a) 



0. (2.21b) 



The common denominator is never zero, and is sufficient to find the zeros in the numerators. 
The solutions 9 = 0,7i satisfy both equations simultaneously giving the solution set A = 
{0,7r} X [0,27r). In the cartesian product the first set gives the possible values of the 9 
coordinate, the second set the possible values of the ip coordinate. 

Setting ip = 0, 7r/2, vr, or 37r/2 solves Eq. (12.21bp : substituting these values into (12. 2 lap 
gives 

sin^cos^(/3 - /i)m^ = 0, = 0, vr, (2.22a) 
sin e cos e (/3 - Is) = 0, ip = 7i/2,37i/2. (2.22b) 



9 



For these equations 9 = 0,7r/2,7r, are solutions. In the same notation used before the new 
solutions are B = {7r/2} x {0,7r}, and C = {71/2} x {7r/2, 37r/2}. 

The complete solution set for the Euler problem is Se = AU B U C. For each of these 
sets is possible to obtain a Vm-m curve simply by evaluating the effective potential on each 
of the solution sets 

2 

Kn(^) = 77^, (2.23a) 

2 

777 

Vm{B) = — , (2.23b) 

All 

2 

771 

VmiC) = — . (2.23c) 

The E-7n diagram consists of three parabolas and the regions enclosed between them. 
This diagram is shown in Fig. [T]for the three molecules considered here (cf. §3.4 of Ref. I^; 



Ch. 14 of Ref. llOl ). The different lines represent the three parabolas (12.231) and correspond 
physically to rotations of the top about the body fixed axes (in a positive or negative sense). 
The region below the red curve is physically inaccessible. In the axis convention used to 
obtain the asymmetric top Hamiltonian the red curve, equation fl2.23bp . corresponds to 
stable rotation about the body fixed x-axis; the green parabola, equation fl2.23cp . to an 
unstable rotation about the y-axis, and the blue curve to a stable rotation about the 2-axis. 

For a given {E,m) point in Fig. [1], the effective potential fl2.17p gives an equation to 
solve in 6 or ip. The solutions are found as contours of constant Vm = E in the O-ijj 
configuration space, the Poisson sphere S"^. These contours define the classically accessible 
(Vm ^ E) and forbidden [Vm > E) regions for given E and m. The topology of the different 
solutions on the (Ojtp) sphere is characterized by the (multi- valued) genus 6|, which for 
every connected classically allowed region counts the number of disjoint discs associated 
with forbidden motion that are removed from the sphere. 

Figure [2] shows the classically accessible regions (in black) of S"^ for iodopentafluoroben- 
zene. Panel |2]|^a) represents an {E,m) point located between the red and green parabolas 
in Fig. [T](c); the multivalued genus is 1,1 since there are two disjoint classically accessible 
regions and from the point of view of each there is one white disc removed from the sphere. 
In the same way, panel [2](b) is the accessible region for a {E,7n) point located between 
the green and blue parabolas of[T](c), but now the genus is 2 since two white discs are re- 
moved. Finally the region above the blue parabola in W^c) has genus since all the sphere 
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is classically accessible (panel [2]^c)). 
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3. E-m DIAGRAM FOR ASYMMETRIC TOP MOLECULES IN STATIC ELEC- 
TRIC FIELDS 



Energy-momentum diagrams for the asymmetric top in an external gravitational field 



have been studied by Katok 



301 1 and by Gashenenko and Richter 5l|. An important con- 



clusion in these works is that, in the study of relative equilibria and their bifurcations in 
asymmetric tops, the relevant parameters are the ratios between two of the moments of 
inertia and the third one, e.g., /1//3 and 12/13- 

In this section we treat as examples three asymmetric top molecules of recent theoretical 



and experimental interest 



13, 



39|: the near-prolate top iodobenzene (CgHsI) 



21 



53j, the 



near-ob 
zene 



ate top pyridazine (C4H4N2^ 



l3, 3; intermediate case iodopentafluoroben- 



21| . Relevant physical parameters for these molecules are given in Table [H 
In Fig. [3] we show the definition of the body fixed frame for these molecules. The moment 
of inertia Jj is related to the rotational constant Bi by Bi = (2Jj)~^. The rotational constants 
Bi, the field parameter u, and the energy are all scaled by B3. For iodobenzene an electric 
field of Es = 25kVcm~^ 53|] gives uj/B^ = 4.52; for pyridazine a field of strength es = 



56kVcm ^ |20| gives u/B-^ = 19.04; in iodopentafluorobenzene, an electric field of £5 
25kVcm-i gives uj/B^ = 18.93 



17 



21 



53|. 



For the asymmetric top molecule in an electric field, the relative equilibria are found by 
solving the equations (setting Auj2 = Au^ = in Equations fl2.18p ): 



+ u sin6 = 0, 



where 



dVm sin 9 cos 9 (Ji cos^ ip + h sin^ ip — Is) 

'W^ [(IC3) ■ C3]2 

dVm sin COS Ip sin^ 9 {I2 ~ 11)171^ 



[\Cs) -03 = 11 sin^ 9 cos^ ip + h sin^ 9 sin^ ip + h cos^ 0. 



(3.1a) 
(3.1b) 

(3.2) 



Again we have the solution set A = {0, vr} x [0, 2n). There are two solution subsets Ai C A, 
and A2 C A, given by Ai = {0} x [0,27r), and A2 = {vr} x [0,27r). The Vm-m curves for 
these subsets are 

i2 



777 ^ 



Co-, 



Vm{A2 



m 
2L 



+ UJ. 



(3.3a) 
(3.3b) 
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The solutions ip = 0, 7r/2, vr, 37r/2 to (IS.lbp must be substituted into fl3.1al) . After rear- 
ranging and dividing by Is, this gives 



for ip = 0,71 and ip = 7r/2,37r/2 respectively, and ii = h/I-s, 12 = h/h- For the case 
Ii > I2 > I3, i.e., ii > 12 > 1, the right hand sides (RHS) of these equations are functions 
with maxima at 6* = 0, vr and a minimum at 6* = 7r/2; the left hand sides (LHS) are functions 
with a minimum at ^ = 0, a maximum at 6 = n, and a fixed zero at ^ = 7r/2. In Fig. HJ the 
LHS and RHS of these equations are plotted in order to show the nature of the solutions. 
As m is increased the amplitude of variation of the solid curve about zero gets larger until 
it intersects the dashed line at 6' = 0, at which point a solution of the equation is obtained. 
Since the curves do not change their shape, there can exist only one solution for each of the 
equations (13. 4p for a given m. This unique intersection gives a value of 6, which together 
with the respective gives the solution of equations (13. 4p . The form of the curves indicates 
that the intersection occurs initially at 6' = and then moves to larger values of 6 as \m\ 
increases, i.e., the solution Ai bifurcates twice, first for the ip = 0,n curve and then for 
the i/j = 7r/2,37r/2 curve. Since the solid line is always zero at = 7r/2, the intersection 
of the two curves cannot go beyond this point, this means that for large values of m the 
asymmetric top in an electric field behaves like an Euler top. 

The bifurcations in the E-m and E-6 diagrams are shown in Fig. [5l For iodobenzene we 
set uj/Bs = 10, for pyridazine u/B^ = 20, and for iodopentafiuorobenzene uj/B^ = 20. 

In the E-m diagrams there are four regions delimited by different color curves. The 
allowed configuration space corresponding to the regions delimited by the red and green 
curves have genus 1,1; those between the green and the blue curves have genus 2; between 
the blue and magenta curves the genus is 1, and above the magenta curve the genus 0. It 
can be seen that the only difference between these diagrams and those for the Euler top is 
the presence of the region with genus 1. 

The curves themselves represent rotations in at constant 9 and ip. Associated 9 values 
for given energy are shown using the curve colors in panels [5l^b), (d), and (f). The red curve 
is a stable rotation with ip = 0,tt and 9; the green curve is unstable with ip = 7r/2, 37r/2 and 
9 in the RHS panels. The blue and magenta curves are associated with degenerate equilibria 




(3.4a) 
(3.4b) 



(1 — ^2) cos 9 



UJI3 [^2 + (1 



(2) cos^ 9] , 



13 



in the coordinate and ^ = (stable, blue) or ^ = vr (unstable, magenta). Physically the 
blue and magenta curves in Fig. [5] represent the situation when the molecule's dipole is 
oriented with the field (blue) and against it (magenta), while at the same time the molecule 
is rotating about its own z-axis with direction and angular velocity given by m. 
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4. E-m DIAGRAMS FOR ASYMMETRIC TOP MOLECULES IN NONRESO- 
NANT LASER FIELDS 



The effective potential for interaction of a molecule with a nonresonant linearly polarized 
laser field (o; = 0) is 

2 

VmiO, i^) = ^^-^ ~ - Au2 sin^ e sin^ - Acu-s cos^ 9. (4.1) 
2(IC3) ■ Ca 

The relative equilibria are obtained by solving 

dVm sin 9 cos 9 (Ji cos^ + h sin^ ip - I3) dV 



89 [(\C,) . 99 

dVm sinip cosip siri^ 9{l2 — Ii)m'^ dV 



0, (4.2a) 



dip [(IC3) ■ Cal^ dip 

where (IC3) ■ C3 is again given by Eq. (13. 2p and 

dV 

'89 
dV 



0, (4.2b) 



— 2Aco'2 sin6'cos6'sin^'?/' + 2Aco'3 sin6'cos^, (4.3a) 



— — = — 2 Au;2 sin ip cos ip sin^ 9. (4.3b) 

dip 

From (I4.2l) - (l4.3p it is seen that 6* = 0, vr are again simultaneous solutions and the first 
solution set is ^ = {0, tt} x [0, 27r). In contrast with the electric field case, the symmetry of 
the laser interaction gives only one Vm-m curve for A, 

2 

77? 

VmiA) = — - Acos. (4.4) 

It is straightforward to see that 6* = 7r/2 is a solution of (I4.2al) . which after substitution into 
(14.2bl) gives 

sin^/^cosV'(/2 - /i)m2 oa • / / n fA ^\ 

r r ■ 2 I 2Aa;2sm7/'cos?/' = 0, (4.5) 

ii cos^ ip + I2 sm ip 

with immediate solutions ip = 0, 7r/2, vr, 37r/2. The second solution set is therefore B = 
{7r/2} X {0, 7r/2, TT, 37r/2}. There are two subsets, Bi C B and B2 C B, given by Bi = 
{it/2} X {0,7r} and B2 = {tt/2} x {7r/2, 37r/2}. For these the Vm-m curves are 

2 

m 

Vm{B,) = — , (4.6a) 

2 

Tfl 

Kn(i32) = — -Aa;2. (4.6b) 
^-'2 
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Now, ip = 0, 71 /2, 71, and 37r/2 are particular solutions of equation (14.2bP : substituting 
these values into equation (14. 2p and dividing by I3 produces 



'1 - ii) = -2ACJ3/3 [ii + (1 - ii) cos^ 
;i - 12) n? = 2{Au2 - Acu3)/3 [i2 + (1 - 



«2) cos^ 6*] ^ 



(4.7a) 
(4.7b) 



for = 0, 71, and ip = 7t/2, 37r/2 respectively, and as before ii = /1//3 and 12 = h/h, 
with ii > i2 > 1. The nature of the solutions depends on the values of Auj2 and Au^. For 
iodobenzene and pyridazine Au;3 > Au!2 > 0, since for these molecules ai < 02 < 03 (cf. 
Table [Tj). This is however not always the case: for pyridine with the same axis convention 



«! < 0:3 < Q!2 



54l |. which implies Auj2 > Au-s > 0. 



For a physically reasonable laser intensity lO^^Wcm ^ applied to iodobenzene we have 
AU2/B3 = 284.6 and Au-^/B^ = 630.5; the same field applied to pyridazine produces 
AUO2IB3 = 225.7 and Auj^/B-i = 228.9; iodobenzene gives AUJ2/B3 = 2280.06 and 
AUJ3/B3 = 4097.94. These dimensionless ratios are much larger than the corresponding 
energy ratios for physically relevant values of the interaction with a static electric field, 
u = 10 and u = 20. As we wish to investigate the interesting dynamical regime in which 
the effects of both fields are of similar magnitude (see next Section), we reduce the intensity 
of the laser by a factor of 10 to obtain the values listed in Table Ull 

The plot of the RHS and LHS of equations (14.71) is shown in Fig. [6l The situation 
is different from the static electric field case. Now the LHS of the equations, the solid 
horizontal line, moves down as m increases. This line intersects the RHS curve (dashed) at 
6 = 0,7r simultaneously for m values 

.1/2 



Ml 
M2 



± 



± 



1-^1 
2(Acj2 - Au3)Is 



1/2 



(4.8a) 
(4.8b) 



l-i2 

7r/2,37r/2, respectively. As the value of m increases the points 



for ip = 0,7T and ip 

of intersection approach the value 6 = 7r/2 symmetrically, i.e., 7r/2 — 6\eft = bright — 7r/2. 
Finally, at m values 



Ml = ±Miii, 

M2 = ±M222, 



(4.9a) 
(4.9b) 
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for ip = 0,7r and ip = 7r/2,37r/2, respectively, the intersection occurs exactly at 6' = 7r/2. 
For larger values of m there is no solution. 

Considering only positive values of m, equations (14. 7p can be solved for cos^ 9 to get 



9 „ m — Ml 
cos^9 = 4.10a 

r, „ m — Mo 

cos^9 = 4.10b 

M2-M2 

for '0 = 0, vr and if) = 37r/2, respectively. The two solutions of each of these equations 
together with the respective ip values give the solution sets C for = 0, vr and T) for 
ip = 7r/2,37r/2. The Vm-Tn curves for these solutions are 

(Mi-Mi)m2 , m-Mi s 

Vm(C) =r - Aujs 4.11a 

^ ^ 2 [/3(m-Mi) + /i(Mi -m)] Mi - Mi ^ ^ 

Vm(V)=—, - AUJ2 + (Alj2 - Aus) 4.11b 

""^ ' 2[h{m-M2) + h{M2-m)\ ^ ' V2 - M2 ^ ^ 

with m G [Ml, Mi] and m G [M2, M2] for the first and second equation respectively. The 
value of Vm{C) at m = Mi is equal to that of Vm{A) at m = Mi, while its value at m = Mi 
is equal to Kn(i3i) at m = Mi; similarly, the value of Vm{T>) at m = M2 equals Vm{A) at 
m = M2, and at m = M2 equals Vm{S2) at m = M2. These results indicate that in the E-m 
diagram the solution set C is a bridge connecting solutions A and Bi. In like fashion the 
solution set V connects the solutions A and 32- 

Finally, the last solution set is obtained from equations (14. 2 p after removing all the 
common factors 



(Ji cos^ + h sin^ ip — I3) 



m 



2 



(/i sin^ 9 cos^ ip + I2 sin^ 9 sin^ tp + I3 cos^ 9^ 



2 



2Au;2 sin^ tp + 2Au-i = 0, (4.12a) 



(I2 — Ii)m? . / , s 

^— 2 - 2AcJ2 = 0. (4.12b) 

(Ji sin^ 9 cos^ ip + I2 sin^ 9 sin^ ip + I'^ cos^ 6*) 

Solutions of these equations for given m are obtained by finding the intersections of the zero 
contours of the LHS of these equations in the 9-ip space. For the case of the laser field, these 
equations can be rearranged to get m? in terms of 9 and ip 



2{Auj2 sin?/'2 - Acjg) 
J3 — Ji cos^ ip — I2 sin^ ip 
2Acu2 
h-h 



{\C.)-C: 



2 



2 



(4.13a) 
(4.13b) 



17 



where the definition (13.20 is used. With the equations written in this form is clear that we 
must find ip such that 

2(Acj2sin2 7/^- Ac^g) _ 2Au2 ^^^^^ 



(4.15) 



(4.16) 



J3 — Ji cos^ ip — h sin tp Ii — I2 
Rearranging and simphfying gives the condition 

1 _ _ ^3 - h 

AiU2 h- 

which can be written in terms of the polarizabihty 

^3 - 0^2 _ h - h 
ai — a2 h — h 

In general this equality is not satisfied for physical parameter values. For example, table [TTl 
lists the LHS and RHS of equation (14.161) for the molecules of interest. 

As condition (14.151) is not fulfilled for the molecules considered here, the last solution set 
is empty, £" = 0. The complete solution for the asymmetric top rotor in a linearly polarized 
laser field is then given byiS^ = AVJ BUCWD. The Vm-m curves for these solutions in the 
case of iodopentafiuorobenzene, iodobenzene and pyridazine are shown in Figs U\ El and [9l 
respectively. In these figures the solution sets are given by different types of curves. The 
solution set A is given by the red curve, Bi by the green curve and B2 by the blue. The 
bridge solution sets C, connecting A with Bi are magenta; the bridge V, connecting A with 
B2, is the cyan curve. 

For iodopentafiuorobenzene and iodobenzene there are nine different regions delimited by 
different curve types with their corresponding genera indicated with arrows in Fig. [HI The 
classically accessible 6-ip configuration space for each of these regions is shown in Fig. [TOl 
The simplest region is characterized by genus 0; in this region the molecule is free to move 
in any possible configuration in 6-ip as can be seen in Fig. [TOT i). There are two regions with 
genus 2: one for high energy and large m, the other for low energy and small m. In the 
high energy region [TOlfa) shows that the molecule is localized in the equatorial region and 
the poles are forbidden. The second region with genus 2 is shown in [TOT c): in this case the 
forbidden regions are on the equator with ip = 0,7r. For genus 4 there is only one region, 
shown in panel fTOT f): the molecule can access neither the poles nor the equatorial regions 
with = 0,7T. In the case of genus 2,2, panel fTOT d). the accessible region is delocalized in ijj 
and restricted in 6^ to a region close, but not on, the poles. Finally, for genus 1,1,1,1 (H^OT f)) 
the motion is highly localized near 9 = ir/A, 37r/4 and ip = 7r/2, 37r/2. 
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The fact that both the inertia and polarizabihty tensors of pyridazine have near-oblate 
symmetry makes it harder to observe the different regions in the E-m diagram, as can be 
seen from Fig. O As /2 ~ I3 and Aco'2 ~ Acus, the Vm-fn curves for the solution sets A and 
B2 are very similar. The magenta bridging solution C connecting A and Bi is easily seen, but 
the cyan V bridge is much harder to see. As for iodopentafluorobenzene and iodobenzene, 
there are 9 different regions which have the same distribution of genus. 
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5. E-m DIAGRAMS FOR ASYMMETRIC TOP MOLECULES IN COLLINEAR 
FIELDS 



We now consider the rotational dynamics of dipolar asymmetric tops in combined static 
(Section [3]) and nonresonant laser (Section Hj) fields. The polarization of the laser field is 
taken to be collinear with the static field, so that m is a conserved quantity. 

A related problem has been analyzed by Tatarinov [g], [sil , who studied the problem of 
the rotation of a rigid body about a fixed point with a potential 

V = PCs-RcM + ^i\C3)-Cs, (5.1) 

where P and p are constants, and Rcm is the location of the center of mass relative to the 
fixed point. In Tatarinov's work E-m diagrams and genera are obtained for the effective 
potential. 

In molecular terms, the first term of the potential (15.11) corresponds to the interaction of 
an electric field along the space-fixed 2;-axis with the dipole moment of the molecule, where 
the dipole moment vector can point in an arbitrary direction in the molecule-fixed frame, 
R, not only along the molecule fixed 2-axis as in the case studied here. The second term 
of (15.11) has no obvious molecular analog, but the dependence on Euler angles is exactly the 
same as for the laser interaction. In fact, the potential energy for an asymmetric top in 
tilted (non-collinear) fields can be written 

V = -n-C3- {AQC,) ■ Cs, (5.2) 

with the definitions f2 = (cu sin /3, 0, u cos (3), and AQ = diag (0, Auj2, Au^). 

The effective potential for the dipolar asymmetric top in collinear fields, where the dipole 
moment points along the molecule- fixed z-axis, is (cf. eq. (I2.17P ) 

2 

Vm{0, ip) = u COS 9 — Auj-2 sin^ 6 siv? ip — Au^ cos^ 9. (5.3) 

2(IC3) ■ C-s 

The relative equilibria are obtained by solving the equations 

dVrr, sin 6 cos 6 {h cos^ ip + h sin^ ip - /g) dV 



89 [(IC3) . CsY 96 



0, (5.4a) 



dVm sinz/' cosz/) sin^ ^(/2 — /i)m^ dV ^ .1 n 
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with 

(ICa) ■ C3 = /i sin^ 9 cos^ + h sin^ 9 sin^ ^ + I3 cos^ ^, (5.5) 

and 

dV 

— — = a; sin 6' — 2Aco'2 sin 9 cos 6^ sin^ + 2Au;3 sin 6' cos 9, (5.6a) 
o9 

dV 

— — = —2Au2 sin cos ^Jj sin^ 6*. (5.6b) 

In both equations (15. 4p the common factor sin 9 gives the simultaneous solution 6' = or 
TT, so the first solution set is ^ = {0,7r} x [0,27r). As for the static electric field case, this 
solution generates two subsets Ai = {0} x [0, 27r) and A2 = {vr} x [0, 27r). The Vm-m curves 
for these solutions are 

2 

771 

Vm{Ai) = — - u - Aus, (5.7a) 

2 

VU^2) = — + iu- Aius. (5.7b) 

As for the laser interaction case, equation fl5.4bl) has solutions ip = 0, 7r/2, vr, and 37r/2, 
which after substitution into (I5.4al) give 

cos 6^(1 -ii)m^ = -13(00 + 2Auj3Cos9) [ii + (1 -ii)cos^6']\ (5.8a) 
cos ^ (1 - is) = -h - 2(Acu2 - Acug) cos ^] [^2 + (1 - ^2) cos^ ^] ^ . (5.8b) 

Figures fTTj [T2] and [T3] show plots of the LHS and RHS of equations (15.81) for iodopentafiu- 
orobenzene {00 /B^ = 20, Acus/^g = 228.006 and Aoo^/Bs = 409.794), iodobenzene 
(lo/B^ = 10, ALO2/B3 = 28.456 and Auo^/B^ = 63.05) and for pyridazine (uj/Bs = 20, 
ALJ2/B3 = 22.56, and Alo^/B^ = 22.88). 

Figures [H] and [121 show the LHS of (15. 8p for iodopentafluorobenzene and iodobenzene 
(respectively) with m = as a dotted line in both cases; it is clear that on each figure that for 
m = there is only one intersection at ^ > 7r/2. As the value of m increases, the dashed curve 
will intersect the solid line at additional points. The RHS of (I5.8al) gives —I3 [uj + 2Au;3) 
at 6' = and —I3 {u — 2Act;3) at ^ = tt; for these two melecules in the fields specified above 
\uj + 2Aco'3| > \uj — 2Aco'3|, which implies that the LHS curve intersects the RHS curve first 
at 9 = 71, and then at 6* = 0. For larger m only the leftmost intersection remains with the 
two intersections with 9 > tx/2 disappearing simultaneously. This situation is repeated for 
the second equation (15.81) . for which the RHS at ^ = 0, tt gives —I3 [oj — 2(Aci;2 — Aa;3)] and 
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— Js [oj + 2(Aco'2 — Acus)] respectively; for the fields employed we have for both molecules 
that \uj -2 (Aa;2 - Ac^s)] > \uj + 2 {5uj2 - Ac^s)]. 

For pyridazine the situation is different. From Fig. [T5](a) it is clear that the RHS curve 
of fl5.8al) is intersected first at 6* = vr, but the two solutions with Q > -n j2 will now disappear 
before the dashed curve intersects the solid curve at ^ = 0. For (15.8bp there is only one 
solution, which moves from 6* = to larger values of Q as the value of m increases. 

Figures [11] to [13] explain the observed bifurcation structure of the solutions of equations 
(15. 8p . The two panels in Figures [TT] and [T2] and panel [T3l^a) show the existence of a nonbi- 
furcating solution and a pair of solutions emerging from a saddle-node bifurcation. A global 
view of the bifurcations is obtained rearranging equations (15. 8p to 

cos 61(1 - ii)m2 + /3(a; + 2Au;3Cos6') [zi + (1 - ii)cos^6']^ = 0, (5.9a) 
cos ^ (1 - %2) + Jg [CU - 2(ACJ2 - Acjg) COS [zg + (1 - ii) cos^ ^] ^ = 0, (5.9b) 

and then plotting the zero contours of the LHS in the Q-m space. These contours are shown 
in Fig. [14] for the molecules treated here. The contours indicate that a solution branch for 
both equations (15. 9p begins at 6' = 0, branching out of the A\ solution. On the other side, 
the right branch from the saddle-node pair begins at ^ = vr, which means that this solution 
emerges from the A2 solution. The left branch of the saddle-node pair is connected only to 
its right partner. 

The last solution set is obtained after removing all common factors in equations (15.40 . 
and rearranging to obtain vn? as function of Q and 

2 2(A^2 sin^ ^ - Ao^g) cos g - a; ^ -3 (v,^c^\ 

m = WTt ~T ~i T ■ 2 ,\ [VC?,) ■ C3] , (5.10a) 

cos U (i3 — ii cos"^ W — h sm ip) 

m' = p^[i\Cs)-Csr, (5.10b) 
-fi — -'2 

with (IC3) ■ C3 given by equation (15. 5p . The solution is obtained by finding all the values of 
9 and ip such that 

2{Au2sm'^ ip — AuJs) cos9 — uj 2Au2 re; m 

— (5-11) 



cos 9 [I3 — Ii cos^ 4' — I2 sin^ -j/^) Ii — I2' 
which can be rearranged and simplified to get 



2 cost' — = - — . (5.12) 

Wi - /2 Aa;2 / Ao;2 ^ ^ 
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In terms of the polarizability this gives 

cos ^ = ^ f ^ - " . (5.13) 

2Au2 \I1-I2 a2-aij ^ ^ 

In contrast to equation (14. 16^ . equation (15. 13^ depends on the fields. For the fields used in 
iodobenzene we obtain 6 = 0.48967r, while for pyridazine there is no solution. 

The range of m for this solution is obtained from equation (I5.10bll . The smallest value 



of m occurs ai ip = n/2, 37r/2, the largest at = 0, vr, giving the values 



h-h 



' 4Acu2 \h-h At^, 



(5.14) 



where i = 1 for = 0, vr, and i = 2 for ip = 71/2,311/2. A solution of (15.101) therefore exists 
for m G [M2, Ml]. Noting equations (15. 8p . and keeping in mind equation (15.111) . we conclude 
that this solution coincides with the solution of ( I5.10ap at m = M2, and with the solution 
of (l5.10bD at m = Mi. 

Using equations (I5.10bp and (I5.14p we obtain 

1/2 



■ 2n ■ 2 I Mi-m (h-h ^ 

sm ^sm V^ = ^— ^(^^^) . (5.15) 



This equation, together with ( I5.12p and ( ]5.10bp can be substituted into the effective potential 
(I2.17P to get the Vm-m curve 

/ 2A \ 

V^ = {m- Ml/2) I J ^1 ] - cos e - Acu3 cos^ ^, (5.16) 



\h-h/ 

where the last two terms are to be evaluated using equation (I5.13p . The first term shows 
a linear dependence on m. In fact, as seen above, this solution is a bridge between the 
solutions of (15. 8p a.t ip = 7r/2,37r/2 and ip = 0,71. Although it has not been considered, for 
negative values of m the solutions are symmetric with respect to reflection on the m = 
line. 

The E-m diagram for iodopentafluorobenzene is shown in Fig. [151 There are 16 different 
regions delimited by the Vm-fn curves. In the figure there are two parabolas corresponding 
to the solution sets Ai (red) and A2 (green), the Ai parabola is always below the A2 as 
should be evident from equations (15. 7p . Emerging from the Ai parabola at m ^ 10 there 
are two curves, one blue the other magenta. These curves correspond to the independent 
branch solutions of equation (15. 9p . the blue curve is the solution for ip = 0,7i, the magenta 
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one for if) = 7r/2,37r/2. These two curves intersect near m ^ 50, and near this intersection 
they are connected by the bridge solution fl5.16p (cyan). Emerging from the A2 curve, near 
m ^ 10, there are also two curves with the same colors as above; these are the two saddle- 
node bifurcation solutions of equation (15. 9p with the blue curve for the = 0,11 solution and 
the magenta curve for the ip = 7r/2, 37r/2 solution. The saddle- node bifurcation are observed 
as cusps near m ^ AO for the blue curve and m ^ 25 for the cyan curve. The complete E-m 
diagram is obtained when negative values of m are considered, the diagram is symmetric 
about the m = axis. For the full diagram the bifurcation partners form "smiles" that are 
typical features of this type of system |3l| (for analogous structures in the case of diatomic 



molecules in coUinear fields, see 



25|). 



The accessible O-ip configuration space for the different regions of Fig. [TSl is shown in Fig. 
[T6l Comparing with the figure for the laser interaction only (Fig. [TOl) . we see that the genus 
2,2, panel [TOT d). is not present when the electric field is turned on; instead, there is a genus 
1,2 region, panel [T6l(l) . This indicates that the effect of the electric field is relatively strong 
compared to the laser and that it tends to align the molecule dipole along the space fixed 
2;- axis. This is also the case for the genus 3 region of Fig. fTGT j). which is absent in Figure 
[TOl Note also that the genus 1,1,1,1 configurations are different in both figures. 

For iodobenzene the E-m diagram. Fig. [TTl is simpler than the one for iodopenfiuoroben- 
zene. Fig. [151 Since the rotational constants for this molecule are considerably larger than 
those for iodopentafluorobenzene, the E-m diagram extends over a smaller energy range. 
There are only few qualitative changes with respect to Fig. [151 The most important differ- 
ence is that the lower smile lies completely within the 2 region and does not go over to the 
1,1 region as for iodopentafluorobenzene. This eliminates the lower 2 region observed in 
the lower panel of Fig. [151 and also the 1,1 region in the same panel. 

The E-m diagrams for pyridazine are shown in Fig. [181 In panel (a), where uj / = 20, 
Auj2/B;i = 22.56, and AuJs/B;^ = 22.88, we obtain a very simple E-m diagram with only 6 
regions. There is only one "smile" with genera 2 and 3. The second smile does not form 
since there is only one branch for the ip = 7r/2,37r/2 solution. In panel (b) we show the 
E-m diagram for the same value of electric field but for a laser field 10 times more intense. 
In this case, the smile has grown bigger, overlapping the regions with genera 1, 2, and 1,1. 
It is observed that below and above the overlapped regions the genus is the same as for 
the weaker laser field case, while inside the smile the genus can vary according to the E-m 
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values. 
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6. SUMMARY AND CONCLUSIONS 



In this paper we have studied the classical mechanics of rotational motion of dipolar 
asymmetric top molecules in strong external fields. Static electric fields, linearly polarized 
nonresonant laser fields, and collinear combinations of the two were investigated. The par- 
ticular asymmetric top molecules iodobenzene, pyridazine and iodopentafluorobenzene have 
been treated for physically relevant field strengths. 

Following Katok 30|] and Tatarinov 311], we have computed diagrams of relative equlibria 



in the E-m plane; the relative equilibria correspond physically to periodic motions with the 
two Euler angles 6 and i/j constant [g], [sqI. We have also examined the classically allowed 
6-ip configuration space for different regions of the E-m diagrams, and have classified the 
configuration space topology according to their genus j^. We anticipate that this classi- 
cal mechanical investigation will be useful in the difficult problem of assigni ng quan tum 



mechanical eigenstates and energy levels for asymmetric tops in external fields 
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TABLE I: Molecular parameters. 



Molecule ,63" Bi/B^ do^ ai"^ 02 as Point group 

CgHgl'^ 0.189 0.117 0.132 1.70 10.2 15.3 21.5 Cg^ 

C4H4N2'= 0.208 0.490 0.970 4.14 5.84 10.29 10.35 C2V 

CeFsl/ 0.034 0.264 0.359 1.54^ 10.5 17.9 23.8 C2„ 



"Rotational contant in cm~^ 
''Dipole moment in Debye 
'^Polarizabilitiasilri"^'^ 
'^Iodobcnzen(F |TrlS3| 
'^Pyridazine [5J, l55 
■'^lo dopentafluorobenzene 
^ab-initio 3-21G 



2l| 
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TABLE II: Hamiltonian parameters and asymmetries with respect to y-axis. 



Molecule u:/B^ Aco2/Bs Au^/Bs ^(^-^-^'^3) 

CfiHsI 10 28.46 63.05 -69.88 -1.22 -6.61 

C4H4N2 20 22.56 22.88 -0.64 -0.014 -0.045 

CeFgl 15 228.0 409.8 -363.58 -0.797 -1.785 
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Figure captions 



FIG. 1: E-m diagrams for free asymmetric top molecules (Euler tops), (a) lodobenzene; (b) 
Pyridazine; (c) lodopentafiuorobenzene. 

FIG. 2: Regions of classically allowed motion in 6-ip configuration space for lodopentafiuo- 
robenzene for the different regions of Figure [T](c). 

FIG. 3: Definition of body-fixed coordinate frames for lodobenzene (upper panel), pyridazine 
(middle panel), and lodopentafiuorobenzene (lower panel). Hydrogen and fluorine atoms are 
represented by white circles, carbon atoms by black circles, and nitrogen and iodine atoms 
by grey circles, respectively. In each case the x-axis is perpendicular to the plane of the 
molecule, and the dipole moment points along the body-flxed 2-axis. 

FIG. 4: Plots of RHS (dashed), and LHS (solid) of equations (13.41) for m values at which 
curves do not intersect. Left panel: equation fl3.4ap : right panel: equation (13.4bp . 
FIG. 5: E-m and 6-m bifurcation diagrams for dipolar asymmetric top molecules in a static 
electric fleld. (a), (b) lodobenzene {uj / = 10); (c), (d) pyridazine {ui/B^ = 20); (e), (f) 
iodopentafluorobenzene {uj / = 20). 

FIG. 6: Plots of RHS (dashed), and LHS (solid) of equations (14.71) for m values at which 
curves do not intersect. Left panel: equation fl4.7al) (zi = 1.25, (1 — ii)m^ = —5, 2AC1J3/3 = 
10); right panel: equation fliJbl) (z2 = 1.2, (1 - ^2)^^ = -4, 2{Au2 - Aljs)^ = 7). 
FIG. 7: E-m diagram for the asymmetric top molecule iodopentafluorobenzene in a nonres- 
onant laser fleld: AUJ2/B3 = 228.006 and Acjg/fig = 409.794. 

FIG. 8: E-m diagram for the asymmetric top molecule lodobenzene in a nonresonant laser 
fleld: AU2/B3 = 28.456 and Au^/B^ = 63.05. 

FIG. 9: E-m diagram for the asymmetric top molecule pyridazine in a nonresonant laser 
field: AUJ2/B3 = 22.56 and AU3/B3 = 22.88. 

FIG. 10: Classically allowed (red) and forbidden (green) 6-ip configuration space and asso- 
ciated genera for the various regions of Figures [7] and [HI 

FIG. 11: RHS (solid line) and LHS for m = 20 (dashed line) of equations (15.81) for 
lodopentafiuorobenzene. (a) Equation (I5.8ap : (b) equation fl5.8bl) . Dotted line: m = 0. 
FIG. 12: RHS (solid line) and LHS for m = 10 (dashed line) of equations (15. 8p for lodoben- 
zene. (a) Equation fl5.8ap : (b) equation fl5.8bp . Dotted line: m = 0. 
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FIG. 13: RHS (solid line) and LHS for m = 4 (dashed line) of equations (15.81) for pyridazine. 

(a) Equation (15. Sal) : (b) equation (]5.8b|) . Dotted line: m = 0. 

FIG. 14: Zero contours of the LHS of equations (15.91) . lodobenzene: (a) eq. (I5.9al) . (b) eq. 
fl5.9bp . Pyridazine: (c) eq. (I5.9al) : (d) eq. fl5.9bp . lodopentafluorobenzene: (e) eq. (I5.9al) . (f) 
eq. IKUh^ . 

FIG. 15: E-m diagram for the dipolar asymmetric top molecule lodopentafluorobenzene in 
collinear static electric and nonresonant laser fields. uj/B^ = 20, Auj2/B^ = 228.006 and 
ALJ3/B3 = 409.7945. 

FIG. 16: Classically allowed (red) and forbidden (green) 6-ip configuration space and asso- 
ciated genera for the regions of Figure [151 Genus not included. 

FIG. 17: E-m diagram for the dipolar asymmetric top molecule iodobenzene in collinear 
static electric and nonresonant laser fields. uj/B^ = 10, AUJ2/B3 = 28.456 and Aus/B^ = 
63.05. 

FIG. 18: E-m diagram for the dipolar asymmetric top molecule pyridazine in collinear static 
electric and nonresonant laser fields, (a) uj/B^ = 20, AUJ2/B3 = 22.56, and AUJ2/B3 = 22.88; 

(b) same electric field but a laser field 10 times more intense. 
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FIGURE 2 (C. A. Arango & G. S. Ezra) 
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FIGURE 3 (C. A. Arango & G. S. Ezra) 
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FIGURE 4 (C. A. Arango & G. S. Ezra) 
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FIGURE 6 (C. A. Arango & G. S. Ezra) 
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FIGURE 7 (C. A. Arango & G. S. Ezra) 
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FIGURE 9 (C. A. Arango & G. S. Ezra) 
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FIGURE 10 (C. A. Arango & G. S. Ezra) 
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FIGURE 11 (C. A. Arango & G. S. Ezra) 
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FIGURE 12 (C. A. Arango & G. S. Ezra) 
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FIGURE 13 (C. A. Arango & G. S. Ezra) 
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FIGURE 14 (C. A. Arango & G. S. Ezra) 
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FIGURE 15 (C. A. Arango & G. S. Ezra) 
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FIGURE 16 (C. A. Arango & G. S. Ezra) 
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FIGURE 17 (C. A. Arango & G. S. Ezra) 
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